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Abstract 

In this paper we study the complexity of factorization of polynomials in the free 
noncommutative ring F(a;i, X2, ■ ■ ■, x„) of polynomials over the field F and noncommuting 
variables xi,X2, ■ ■ ■ ,Xn- Our main results are the following: 

• Although V{xi,... ,Xn) is not a unique factorization ring, we note that variahle- 
disjoint factorization in F(a:i,..., Xn) has the uniqueness property. Furthermore, we 
prove that computing the variable-disjoint factorization is polynomial-time equiva¬ 
lent to Polynomial Identity Testing (both when the input polynomial is given by an 
arithmetic circuit or an algebraic branching program). We also show that variable- 
disjoint factorization in the black-box setting can be efficiently computed (where the 
factors computed will be also given by black-boxes, analogous to the work [KT 90 ) 
in the commutative setting). 

• As a consequence of the previous result we show that homogeneous noncommuta¬ 
tive polynomials and multilinear noncommutative polynomials have unique factor¬ 
izations in the usual sense, which can be efficiently computed. 

• Finally, we discuss a polynomial decomposition problem in F(a;i,... ,Xn) which is 
a natural generalization of homogeneous polynomial factorization and prove some 
complexity bounds for it. 


1 Introduction 

Let F be any field and X = {xi, 0:2, • • •, x„} be a set of n free noncommuting variables. 
Let X* denote the set of all free words (which are monomials) over the alphabet X with 
concatenation of words as the monoid operation and the empty word e as identity element. 

The free noncommutative ring F {X) consists of all finite F-linear combinations of mono¬ 
mials in X*, where the ring addition -|- is coefficient-wise addition and the ring multiplication 
* is the usual convolution product. More precisely, let /, g G F(A) and let f{m) G F denote 
the coefficient of monomial m in polynomial /. Then we can write / = Ylm and 

g = ™ product polynomial fg for each monomial m we have 

faim) = X 

m\m 2 =m 
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The degree of a monomial m € X* is the length of the monomial m, and the degree 
deg / of a polynomial / G IF {X) is the degree of a largest degree monomial in / with nonzero 
coefficient. For polynomials f,g € F(X) we clearly have deg(/ 5 ') = deg / + deg s'. 

A nontrivial factorization of a polynomial / G F {X) is an expression of / as a product 
f = gh of polynomials g,h ^ F(X) such that degs > 0 and degh > 0. A polynomial 
/ G F(X) is irreducible if it has no nontrivial factorization and is reducible otherwise. For 
instance, all degree 1 polynomials in F {X) are irreducible. Clearly, by repeated factorization 
every polynomial in F(X) can be expressed as a product of irreducibles. 

In this paper we study the algorithmic complexity of polynomial factorization in the free 
ring F(X). 

Polynomial Factorization Problem 

The problem of polynomial factorization in the commutative polynomial ring F[xi, X 2 ,..., 
is a classical well-studied problem in algorithmic complexity culminating in Kaltofen’s cele¬ 
brated efficient factorization algorithm |Ka89] . Kaltofen’s algorithm builds on efficient algo¬ 
rithms for univariate polynomial factorization; there are deterministic polynomial-time algo¬ 
rithms over rationals and over fields of unary characteristic and randomized polynomial-time 
over large characteristic fields (the textbook [GG] contains a comprehensive excellent treat¬ 
ment of the subject). The basic idea in Kaltofen’s algorithm is essentially a randomized reduc¬ 
tion from multivariate factorization to univariate factorization using Hilbert’s irreducibility 
theorem. Thus, we can say that Kaltofen’s algorithm uses randomization in two ways: the 
first is in the application of Hilbert’s irreducibility theorem, and the second is in dealing 
with univariate polynomial factorization over fields of large characteristic. In a recent paper 
Kopparty et al [KSS14] have shown that the first of these requirements of randomness can 
be eliminated, assuming an efficient algorithm as subroutine for the problem of polynomial 
identity testing for small degree polynomials given by circuits. More precisely, it is shown in 
[KSSI4] that over finite fields of unary characteristic (or over rationals) polynomial identity 
testing is deterministic polynomial-time equivalent to multivariate polynomial factorization. 

Thus, in the commutative setting it turns out that the complexity of multivariate poly¬ 
nomial factorization is closely related to polynomial identity testing (whose deterministic 
complexity is known to be related to proving superpolynomial size arithmetic circuit lower 
bounds). 

Noncommutative Polynomial Factorization 

The study of noncommutative arithmetic computation was initiated by Nisan [N91j in which 
he showed exponential size lower bounds for algebraic branching programs that compute the 
noncommutative permanent or the noncommutative determinant. Noncommutative polyno¬ 
mial identity testing was studied in [BW051 IRS05| . In |BW05| a randomized polynomial 
time algorithm is shown for identity testing of polynomial degree noncommutative arith¬ 
metic circuits. For algebraic branching programs |RS05| give a deterministic polynomial-time 
algorithm. Proving superpolynomial size lower bounds for noncommutative arithmetic cir¬ 
cuits computing the noncommutative permanent is open. Likewise, obtaining a deterministic 
polynomial-time identity test for polynomial degree noncommutative circuits is open. 

In this context, it is interesting to ask if we can relate the complexity of noncommutative 
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factorization to noncommutative polynomial identity testing. However, there are various 
mathematical issues that arise in the study of noncommutative polynomial factorization. 

Unlike in the commutative setting, the noncommutative polynomial ring F(X) is not a 
unique factorization ring. A well-known example is the polynomial 

xyx + X 

which has two factorizations: x{yx + 1) and {xy + l)x. Both xy + 1 and yx + 1 are irreducible 
polynomials in F(A). 

There is a detailed theory of factorization in noncommutative rings [Co85t I Co] . We will 
mention an interesting result on the structure of polynomial factorizations in the ring R = 
F(A). 

Two elements a,b & R are similar if there are elements a', b' £ R such that aV = a'b, and 
(i) a and a' do not have common nontrivial left factors, (ii) b and b' do not have common 
nontrivial right factors. 

Among other results, Cohn |Co] has shown the following interesting theorem about fac¬ 
torizations in the ring R = F(A). 

Theorem 1.1 (Cohn’s theorem) For a polynomial a G F(A) let 

a = 0102 ... Or- and a = 6162 ■ ■ - bg 

be any two factorizations of a into irreducible polynomials in F(X). Then r = s, and there is 
a permutation vr of the indiees { 1 , 2 ,..., r} such that a-i and are similar polynomials for 
1 < i < r. 

For instance, consider the two factorizations of xyx + x above. We note that polynomials 
xy+1 and yx + 1 are similar. It is easy to construct examples of degree d polynomials in 
F(X) that have distinct factorizations. Cohn [Co85| discusses a number of interesting 
properties of factorizations in F(A). But it is not clear how to algorithmically exploit these 
to obtain an efficient algorithm in the general case. 

Our Results 

In this paper, we study some restricted cases of polynomial factorization in the ring F(X) 
and prove the following results. 

• We consider variable-disjoint factorization of polynomials in F(A) into variable-disjoint 
irreducibles. It turns out that such factorizations are unique and computing them is 
polynomial-time equivalent to polynomial identity testing (for both noncommutative 
arithmetic circuits and algebraic branching programs). 

• It turns out that we can apply the algorithm for variable-disjoint factorization to two 
special cases of factorization in F(A): homogeneous polynomials and multilinear poly¬ 
nomials. These polynomials do have unique factorizations and we obtain efficient algo¬ 
rithms for computing them. 

• We also study a natural polynomial decomposition problem for noncommutative poly¬ 
nomials and obtain complexity results. 
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2 Variable-disjoint Factorization Problem 

In this section we consider the problem of factorizing a noncommutative polynomial / € F(X) 
into variable disjoint factors. 

For a polynomial / € F(X) let Var{f) C X denote the set of all variables occurring in 
nonzero monomials of /. 

Definition 2.1 A nontrivial variable-disjoint factorization of a polynomial f € F(X) is a 
factorization 

f = gh 

such that degg > 0 and degh > 0, and Var{g) n Var{h) = 0. 

A polynomial f is variable-disjoint irreducible if does not have a nontrivial variable- 
disjoint factorization. 

Clearly, all irreducible polynomials are also variable-disjoint irreducible. But the converse 
is not true. For instance, the familiar polynomial xyx -)- x is variable-disjoint irreducible 
but not irreducible. Furthermore, all univariate polynomials in F(X) are variable-disjoint 
irreducible. 

We will study the complexity of variable-disjoint factorization for noncommutative polyno¬ 
mials. First of all, it is interesting that although we do not have the usual unique factorization 
in the ring F(X), we can prove that every polynomial in F(X) has a unique variable-disjoint 
factorization into variable-disjoint irreducible polynomialso 

We can exploit the properties we use to show uniqueness of variable-disjoint factorization 
for computing the variable-disjoint factorization. Given / G F(X) as input by a noncommu¬ 
tative arithmetic circuit the problem of computing arithmetic circuits for the variable-disjoint 
irreducible factors of / is polynomial-time reducible to PIT for noncommutative arithmetic 
circuits. An analogous result holds for / given by an algebraic branching programs (ABPs). 
Hence, there is a deterministic polynomial-time algorithm for computing the variable-disjoint 
factorization of / given by an ABP. Also in the case when the polynomial / € F(A) is 
given as input by a black-box (appropriately defined) we give an efficient algorithm that gives 
black-box access to each variable-disjoint irreducible factor of /. 

Remark 2.2 Factorization of commutative polynomials into variable-disjoint factors is stud¬ 
ied by Shpilka and Volkovich in [SVlOj . They show a deterministic reduction to polynomial 
identity testing. However, the techniques used in their paper are specific to commutative 
rings, involving scalar substitutions, and do not appear useful in the noncommutative case. 
Our techniques for factorization are simple, essentially based on computing left and right 
partial derivatives of noncommutative polynomials given by circuits or branching programs. 

2.1 Uniqueness of variable-disjoint factorization 

Although the ring F(A) is not a unique factorization domain we show that factorization into 
variable-disjoint irreducible factors is unique. 

For a polynomial / € F(X) let mon{f) denote the set of nonzero monomials occurring in 

/. 

^Uniqueness of the factors is upto scalar multiplication. 
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Lemma 2.3 Let f = gh such that Var{g) D Var{h) = (p and \Var{g)\, \Var{h)\ > 1. Then 

mon{f) = {mw\m G mon{g),w G mon{h)}. 

Moreover, the coefficient of mw in f is the product of the coefficients of m in g and w in h. 

Proof. Let m G mon{g) and w G mon{h). We will argue that the monomial mw is in mon{f) 
and can be obtained in a unique way in the product gh. Namely, by multiplying m G mon{g) 
with w G mon{h). 

Assume to the contrary that for some u G mon{g) and v G mon{h), u ^ m we have 
mw = uv. Note that m and u are words over Var{g) and w,v are words over Var{h). 
Clearly, |m| 7 ^ |u| because m ^ u. Without loss of generality, we assume that |u| > \m\. As 
mw = uv, it follows that for some word s G X*, |s| > 0, we have u = ms and w = sv. 

That implies Var{s) C Var{g) n Var{h) which contradicts variable disjointness of g and 
h. □ 

Lemma 2.4 Let f = g.h and f = u.v be two nontrivial variable-disjoint factorizations of f. 
That is, 


Var{g) fl Var{h) = (p 
Var{u) n Var{v) = (p. 

Then either Var{g) C Var{u) and Var{h) D Var{v) or Var{u) C Var{g) and Var{v) A 
Var{h). 

Proof. Suppose to the contrary that x G Var{g) \ Var{u) and y G Var[u) \ Var{g). Let m 
be any monomial of / in which variable y occurs, and let m = mim 2 such that mi G mon{g) 
and m 2 G mon{h). Clearly, y must occur in m 2 . Similarly, for any monomial m' of / in 
which variable x occurs, if m' = m'^m^, for m'l G mon{g) and m^ G mon(h), then x must 
occur in m[. Thus, x G Var{g) and y G Var{h). Similarly, it also holds that y G Var{u) and 
X G Var{v). Thus we have: 

1. X G Var{g) and y G Var{h) 

2. y G Var{u) and x G Var{v). 

Hence, there are monomials of / in which both x and y occur. Furthermore, (1) above 
implies that in each monomial of / containing both x and y, x always occurs before y. On 
the other hand, (2) implies y occurs before x in each such monomial of /. This contradicts 
our assumption. □ 

Lemma 2.5 Let f G F(A) and suppose f = gh and f = uv are two variable-disjoint factor¬ 
izations of f such that Var{g) = Var{u). Then g = au and h = fdv for scalars a, /3 G F. 

Proof. As f = gh = uv are variable-disjoint factorizations, by Lemma 12.31 each monomial 
m G mon{f) is uniquely expressible as a product m = mim 2 for mi G mon{g) and m 2 G 
mon{h), and as a product m = m'j^m^ for m'^ G mon{u) and m^ G mon(v). 
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As Var{g) = Var{u), we notice that Var{h) = Var{v), because Var{f) = Var{g) l+l 
Var{h) and Var{f) = Var{u) l±) Var{v). Now, from m = mim 2 = m'^m '2 it immediately 
follows that mi = m'l and m 2 = m^. Hence, mon{g) = mon{u) and mon{h) = mon{v). 

Furthermore, if m is a monomial of maximum degree in g, by taking the left partial 
derivative of / w.r.t. m we have 


dm 


a'h = P'v, 


where a' and j5' are coefficient of m in g and u respectively. It follows that h = j5v for 
some /? G F. Similarly, by taking the right partial derivative of / w.r.t. a maximum degree 
monomial in h we can see that g = au for a G F. □ 


We now prove the uniqueness of variable-disjoint factorizations in F(A). 


Theorem 2.6 Every polynomial in F(Ai) has a unique variable-disjoint factorization as a 
product of variable-disjoint irreducible factors, where the uniqueness is upto scalar multiples 
of the irreducible factors. 


Proof. We prove the theorem by induction on the degree d of polynomials in F(W). The 
base case d = 1 is obvious, since degree 1 polynomials are irreducible and hence also variable- 
disjoint irreducible. Assume as induction hypothesis that the theorem holds for polynomials 
of degree less than d. Let / G F(X) be a polynomial of degree d. If / is variable-disjoint 
irreducible there is nothing to prove. Suppose / has nontrivial variable-disjoint factors. Let 
f = gh he a nontrivial variable-disjoint factorization such that Var(g) is minimum. Such 
a factorization must exist because of Lemmas 12.31 and 12.41 Furthermore, the set Var[g) is 
uniquely defined and by Lemma [23] the left factor g is unique upto scalar multiples. Applying 
induction hypothesis to the polynomial h now completes the induction step. □ 


2.2 Equivalence with PIT 

Theorem 2.7 Let f G F(Ai) be a polynomial as input instance for variable-disjoint factor¬ 
ization. Then 

1. If f is input by an arithmetic circuit of degree d and size s there is a randomized 
poly(s, d) time algorithm that factorizes f into variable-disjoint irreducible factors. 

2. If f is input by an algebraic branching program there is a deterministic polynomial-time 
algorithm that factorizes f into its variable-disjoint irreducible factors. 

Proof. We first consider the most general case of the algorithm when / G F(A) is given by an 
arithmetic circuit of polynomially bounded syntactic degree. The algorithm specializes to the 
other cases too. The algorithm must compute an arithmetic circuit for each variable-disjoint 
factor of /. We explain the polynomial-time reduction to Polynomial Identity Testing (PIT) 
for noncommutative arithmetic circuits. 

Let d = deg/ which can be computed by first computing circuits for the polynomially 
many homogeneous parts of /, and then using PIT as subroutine on each of them. 

Next, we compute a monomial m G mon{f) of maximum degree d with polynomially many 
subroutine calls to PIT. The basic idea here is to do a prefix search for m. More precisely, 
suppose we have computed a prefix m' such that the left partial derivative computes a 
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polynomial of degree d — \m'\ then we extend m' with the first variable Xi ^ X such that 
computes a polynomial of degree d — \m'\ — 1. Proceeding thus, we will compute the 
lexicographically first monomial m € mon{f) of degree d. 

Now, starting with di = 1 we look at all factorizations m = mi.m 2 of monomial m with 
|mi| = di. 

Let h = and g = g;^, circuits for which can be efficiently computed from the given 
circuit for /. Let a, fd and 7 be the coefficients of m in /, mi in g, and m 2 in h, respectively 
(which can be computed in deterministic polynomial time |AMS10] 1. 

Next we compute Var{g) and Var{h). Notice that xi 0 Var{g) if and only if g and g\^.^Q 
are not identical polynomials. Thus, we can easily determine Var{g) and Var{h) with n 
subroutine calls to PIT. 

Clearly, if Var{g) n Var{h) ^ 0 then we do not have a candidate variable-disjoint factor¬ 
ization that splits m as mim 2 and we continue with incrementing the value of di. Else, we 
check if 


a 


f = ^gh, 

Pi 


( 1 ) 


with a subroutine call to PIT. If / = -^gh then ^g is the unique leftmost variable- 
disjoint irreducible factor of /(upto scalar multiplication), and we continue the computation. 
Otherwise, we continue the search with incrementing the value of di. 

In the general step, suppose we have already computed / = gig 2 ■ ■ ■ gihi, where gi, g 2 , ■ ■ ■, gi 
are the successive variable-disjoint irreducible factors from the left. There will be a corre¬ 
sponding factorization of the monomial m as 


m = mim 2 


. mimi 


where mj occurs in 5 ^, 1 < j < i and m[ occurs in hi. Notice that the polynomial 

~ dmimi m- Computed by a small noncommutative arithmetic circuit obtained 

from / by partial derivatives. This describes the overall algorithm proving that variable- 
disjoint factorization of / is deterministic polynomial-time reducible to PIT when / is given 
by arithmetic circuits. The algorithm outputs arithmetic circuits for each of the variable- 
disjoint irreducible factors of /. 

We note that the algorithm specializes to the case when / is given by an ABP. The 
variable-disjoint irreducible factors are computed by ABPs in this case. 

Finally, to complete the proof we note that in the case of noncommutative arithmetic cir¬ 
cuits of polynomial degree there is a randomized polynomial time PIT algorithm [BW05| . For 
polynomials given as ABPs, there is a deterministic polynomial-time PIT algorithm [RS05] . 

□ 


Next we consider variable-disjoint factorization of polynomials input in sparse representa¬ 
tion and show that the problem is solvable in deterministic logspace (even by constant-depth 
circuits). Recall that AC*^ circuits mean a family of circuits {C„}, where for length 
n inputs, such that Cn has polynomially bounded size and constant-depth and is allowed 
unbounded fanin AND and OR gates. The class of TC*^ circuits is similarly defined, but 
is additionally allowed unbounded fanin majority gates. The logspace uniformity condition 
means that there is a logspace transducer that outputs Cn on input 1 ” for each n. 


Theorem 2.8 Let f € F(X) be a polynomial input instance for variable-disjoint faetorization 
given in sparse representation. 
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(a) When ¥ is a fixed finite field the variable-disjoint factorization is computable in deter¬ 
ministic logspace (more precisely, even by logspace-uniform AC^ circuits). 

(b) When F is the field of nationals the variable-disjoint factorization is computable in de¬ 
terministic logspace (even by logspace-uniform TC^ circuits). 

Proof. We briefly sketch the simple circuit constructions. Let us first consider the case when 
F is a fixed finite field. The idea behind the AC^ circuit for it is to try all pairs of indices 
l<a< 6 <din parallel and check if / has an irreducible factor of degree 6 — a + 1 located in 
the (a, b) range. More precisely, we will check in parallel for a variable-disjoint factorization 
/ = 919293 , where <71 is of degree a — 1 , 5^2 is variable-disjoint irreducible of degree b — a + 1, 
and gs is of degree d — b. If such a factorization exists then we will compute the unique (upto 
scalar multiples) irreducible polynomial 53 located in the (a, b) range of degrees. 

In order to carry out this computation, we take a highest degree nonzero monomial m of 
/ and factor it as m = mim 2 m 3 , where m 2 is located in the (a, 6 ) range of m. Clearly, the 
polynomial 



has to be <72 upto a scalar multiple. Likewise, is 93 and g^J^^ is 91 upto scalar 

multiples. Since / is given in sparse representation, we can compute these partial derivatives 
for each monomial in parallel in constant depth. 

We can then check using a polynomial size constant-depth circuit if / = 919293 upto an 
overall scalar multiple (since we have assumed F is a constant size finite field). For a given a, 
the least b for which such a factorization occurs will give us the variable-disjoint irreducible 
factor located in the (a, b) range. For each a < d we can carry out this logspace computation 
looking for the variable-disjoint irreducible factor (if there is one) located at (a, b) range for 
some b > a. 

In this manner, the entire constant-depth circuit computation will output all variable- 
disjoint irreducible factors of / from left to right. This concludes the proof for the constant 
size finite field case. 

Over the field of rationals, we use the same approach as above. The problem of checking 
if / = 919293 upto an overall scalar multiple will require integer multiplication of a constant 
number of integers. This can be carried out by a constant-depth threshold circuit. This gives 
the TC^ upper bound. □ 

We now observe that PIT for noncommutative arithmetic circuits is also deterministic 
polynomial-time reducible to variable-disjoint factorization, making the problems polynomial¬ 
time equivalent. 

Theorem 2.9 Polynomial identity testing for noncommutative polynomials f € F(X) given 
by arithmetic circuits (of polynomial degree) is deterministic polynomial-time equivalent to 
variable-disjoint factorization of polynomials given by noncommutative arithmetic circuits. 

Proof. Clearly proof of Theorem 12.71 gives a reduction from variable-disjoint factorization 
to PIT. To see the other direction, let / G F(A) and y,z 0 A be two new variables. We 
observe that the polynomial / G F(A) is identically zero iff the polynomial / -|- y.z has a 
nontrivial variable-disjoint factorization. This gives a reduction from PIT to variable-disjoint 
factorization. □ 






2.3 Black-Box Variable-Disjoint Factorization Algorithm 

In this subsection we give an algorithm for variable-disjoint factorization when the input 
polynomial / G F {X) is given by black-box access. We explain the black-box model below: 

In this model, the polynomial / G F(A) can be evaluated on any n-tuple of matrices 
(Ml, M2,..., Mn) where each Mi \s a. t x t matrix over F (or a suitable extension field of F) 
and get the resulting t x t matrix /(Mi, M2,..., M„) as output. 

The algorithm for black-box variable-disjoint factorization takes as input such a black-box 
access for / and outputs black-boxes for each variable-disjoint irreducible factor of /. More 
precisely, the factorization algorithm, on input i, makes calls to the black-box for / and works 
as black-box for the variable-disjoint irreducible factor of /, for each i. 

The efficiency of the algorithm is measured in terms of the number of calls to the black-box 
and the size t of the matrices. In this section we will design a variable disjoint factorization 
algorithm that makes polynomially many black-box queries to / on matrices of polynomially 
bounded dimension. We state the theorem formally and then prove it in the rest of this 
section. 

Theorem 2.10 Suppose f G F(A) is a polynomial of degree bounded by D, given as input 
via black-box access. Let f = / 1/2 ■ ■ ■ fr be the variable-disjoint factorization of f. Then there 
is a polynomial-time algorithm that, on input i, computes black-box access to the i*'^ factor fi. 

Proof. In the proof we sketch the algorithm to give an overview. We then explain the 
main technical parts in the next two lemmas. The black-box algorithm closely follows the 
white-box algorithm already described in the first part of this section. The main steps of the 
algorithm are: 

1. We compute a nonzero monomial m of highest degree d < D in f by making black-box 
queries to /. We can also compute the coefficient f{m) of m by making black-box 
queries to /. 

2 . Then we compute the factorization of this monomial 

m = 77117712 • • • rUr 

which, as in the white-box case, is precisely how the monomial m factorizes according 
to the unique variable-disjoint factorization / = /1/2 ... /^ of / (that is nii G mon{fi)). 

3 . In order to find this factorization of the monomial m, it suffices to compute all factor¬ 
izations 777 = m'm" of 777 such that / = (^^)(^^) is a variable disjoint factorization 
(this is done using Equation [U exactly as in the proof of Theorem 12 . 7 p . 

Finding all such factorizations will allow us to compute the 777*. However, since we have 
only black-box access to / we will achieve this by creating and working with black-box access 
to the partial derivatives and for each factorization 777 = m'm" . We explain the 

entire procedure in detail. 

Starting with = 1 we look at all factorizations m = m'm" of monomial m with \m'\ = di. 

Suppose we have computed black-boxes for h = and g = Let a, ft and 7 be the 
coefficients of m in /, m' in g, and m" in h, which can be computed using black-box queries 
to /, g and h, respectively. 


9 


Next we compute Var{g) and Var{h). Notice that Xi 0 Var{g) if and only if g and g\^.^Q 
are not identical polynomials. This amounts to performing PIT with black-box access to g. 
Thus, we can easily determine Var{g) and Var{h) by appropriately using black-box PIT. 

Clearly, if Var{g) n Var{h) ^ 0 then we do not have a candidate variable-disjoint factor¬ 
ization that splits m as mim 2 and we continue with incrementing the value of di. Else, we 
check if 


/ = ( 2 ) 

with a subroutine call to PIT (like Equation [1] in Theorem 12.7p . li f = -^gh then is 
the unique leftmost variable-disjoint irreducible factor of / (upto scalar multiplication), and 
we continue the computation. Otherwise, we continue the search with incrementing the value 
of di- 

Suppose we have computed the monomial factorization m = mim 2 ■ ■ ■ Triim'-. There is a 
corresponding factorization of / as / = / 1/2 ... fihi, where fi, f 2 , ■ ■ ■, fi are the successive 
variable-disjoint irreducible factors from the left, and rrij occurs in fj, 1 < j < i and m[ 
occurs in hi. Then we have black-box access to hi = ^———: and we can find its leftmost 
irreducible variable-disjoint factor as explained above. This completes the overall algorithm 
sketch for efficiently computing black-boxes for the irreducible variable-disjoint factors of / 
given by black-box. □ 


Remark 2.11 The variable-disjoint irreducible factors of f are unique only upto scalar mul¬ 
tiples. However, we note that the algorithm in Theorem \2.10\ computes as black-box some 
fixed scalar multiple for each variable-disjoint irreducible factor fi. Since we can compute the 
coefficient of the monomial m in f (where m is computed in the proof of Theorem \2.10\) . we 
can even ensure that the product / 1/2 ... /r equals f by appropriate scaling. 


Lemma 2.12 Given a polynomial f of degree d with black-box access, we can compute a 
degree-d nonzero monomial of f, if it exists, with at most nd calls to the black-box on {D + 
l)2d X {D -\- l)2d matrices. 


Proof. Similar to the white-box case, we will do a prefix search for a nonzero degree-d 
monomial m. We explain the matrix-valued queries we will make to the black-box for / in 
two stages. In the first stage the query matrices will have the noncommuting variables from 
X as entries. In the second stage, which we will actually use to query /, we will substitute 
matrices (with scalar entries) for the variables occurring in the query matrices. 

In order to check if there is a non-zero degree d monomial in / with Xi as first variable, 
we evaluate f on {D + 1) x {D + 1) matrices Xi,..., Xn defined as follows. 




/ 0 1 0 
: 0 xi 


0 \ 

0 


: ■ • ■ • 0 xi 

VO. 0 y 
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and for i > 1 


/ 0 0 . 0 \ 

: 0 Xi ■ ■. : 

0 

: ■ • • ■ •. 0 Xi 

VO. 0; 


Therefore, f{Xi,... ,Xn) is a (D + l) x (D + l) matrix with entries in F(X). In particular, 
by our choice of the matrices above, the {l,d + 1 )*^ entry of f{Xi,... ,Xn) will precisely be 
where fd is the degree d homogeneous part of /. To check if is nonzero, we can 
substitute random txt matrices for each variable Xi occurring in the matrices Xi and evaluate 
/ to obtain a {{D + l)t) x {{D + l)t) matrix. In this matrix {l,d + 1)*^ block of size txt 
will be nonzero with high probability for t = 2d [BWn5| . In general, suppose m' is a degree 
d monomial that is nonzero in /. We search for the first variable Xi such that m'xi is a 
prefix of some degree d nonzero monomial by similarly creating black-box access to 
and then substituting random txt matrices for the variables and evaluating / to obtain 
a {{D + l)t) X {{D + l)t) matrix. The {l,d + 1 )*^ block of size t x t is nonzero with high 
probability [BW05| if m'xi is a prefix for some nonzero degree d monomial. Continuing thus, 
we can obtain a nonzero monomial of highest degree d in /. □ 

We now describe an efficient algorithm that creates black box access for left and right 
partial derivatives of / w.r.t. monomials. Let m € X* be a monomial. We recall that the left 
partial derivative of / w.r.t. m is the polynomial 


dm 




Similarly, the right partial derivative of / w.r.t. m is the polynomial 


dm 


yy f{m'm)m'. 


Lemma 2.13 Given a polynomial f of degree d with black-box access, there are efficient 
algorithms that give black-box access to the polynomials for any monomials mi, m 2 G 

X*. Furthermore, there is also an efficient algorithm giving black-box access to the polynomial 

(kl\ 

dmi ydm2 J ' 

Proof. Given the polynomial / and a monomial mi, we first show how to create black-box 
access to Let the monomial mi be ... Xi^.. Define matrices Xi ,..., X^ as follows: 
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For 1 <r <k, the entry Tj[r] = 1 if Xi^ = Xi and Tj[r] = 0 otherwise. 

f 

The black-box access to on input matrices Mi,..., of size t x t can be created 

as follows. Note that in the {D + 1) x {D + 1) matrix f{Xi,... ,Xn), the (l,j -|- 1)*^ lo¬ 
cation contains where fj is the degree homogeneous part of /. Now, suppose for 
each variable Xi in the matrix Xi we substitute a t x t scalar-entry matrix Mi and compute 
the resulting /(Xi,... ,Xn) which is now a (D -|- 1 ) x (D -|- 1) block matrix whose entries 


are t xt scalar-entry matrices. Then the block matrix located at the (l,j -|- 1)^* entry for 

f ■ 

jG{2,...,(i-|-l}is the evaluation of on Mi,..., Mn- We output the sum of these matrix 
entries over 2 < j < d -|- 1 as the black-box output. 



The matrix X^^'^ is a (j -|- 1) x {j -|- 1) matrix where, in the figure above the top left block 
is of dimension (j — k) x {j — k + 1). Here Ti[r] = 1 if Xi^ = Xi and Tj[r] = 0 otherwise. 

Finally, define the block diagonal matrix Xi = diag{x\^\ ..., 

We now describe the black-box access to Let Xi ^ Mi be the input assignment of 

t X t matrices Mi to the variables Xi. This results in matrices Xi,l < i < n, where Xi is 
obtained from Xi by replacing Xi by Mi. 

The query /(Xi, X 2 ,..., X„) to / gives a block diagonal matrix diag{Ni^, X^+i,..., Nd). 
Here the matrix Nj is a (j -|- 1) x (j -|- 1) block matrix with entries that are t x t matrices 
(over F). The (1, j -|- 1)*^ block in Nj is evaluated at Mi, M 2 ,..., Mn- 
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Hence the sum of the (1,j + 1)*^ blocks of the different Nj,k < k < D gives the desired 
black-box access to 

□ 

3 Factorization of Multilinear and Homogeneous Polynomials 

In this section we briefly discuss two interesting special cases of the standard factorization 
problem for polynomials in F(X). Namely, the factorization of multilinear polynomials and 
the factorization of homogeneous polynomials. It turns out, as we show, that factorization 
of multilinear polynomials coincides with their variable-disjoint factorization. In the case of 
homogeneous polynomials, it turns out that by renaming variables we can reduce the problem 
to variable-disjoint factorization. 

In summary, multilinear polynomials as well as homogeneous polynomials have unique 
factorizations in F(X), and by the results of the previous section these can be efficiently 
computed. 

A polynomial / G F {X) is multilinear if in every nonzero monomial of / every variable in X 
occurs at most once. We begin by observing some properties of multilinear polynomials. Let 
Var(/) denote the set of all indeterminates from X which appear in some nonzero monomial 
of /. 

It turns out that factorization and variable-disjoint factorization of multilinear polynomials 
coincide. 

Lemma 3.1 Let f € F(A) be a multilinear polynomial and f = gh be any nontrivial factor¬ 
ization of f. Then, Var{g) PI Var(h) = 0. 

Proof. Suppose Xi € Var{g) fl Var{h) for some xi € X. Let mi be a monomial in g of 
maximal degree which also has the indeterminate Xi occurring in it. Similarly, let monomial 
m 2 be of maximal degree in h with Xi occurring in it. The product monomial mim 2 is not 
multilinear and it cannot be nonzero in /. This monomial must therefore be cancelled in 
the product gh, which means there are nonzero monomials of g and m^ of h such that 
mim 2 = m[m 2 . Since deg(m'^) < deg(mi) and deg(m 2 ) < deg(m 2 ) the only possibility is that 
mi = m'l and m 2 = m^ which means the product monomial mim 2 has a nonzero coefficient 
in /, contradicting the multilinearity of /. This completes the proof. □ 

Thus, by Theorem l2.6l multilinear polynomials in F(A1) have unique factorization. Further¬ 
more, the algorithms described in Section [2]2] can be applied to efficiently factorize multilinear 
polynomials. 

We now briefly consider factorization of homogeneous polynomials in F(A1). 

Definition 3.2 A polynomial f € F(A) is said to be homogeneous of degree d if every 
nonzero monomial of f is of degree d. 

Homogeneous polynomials do have the unique factorization property. This is attributed to 
J.H. Davenport in |Ca m- However, we argue this by reducing the problem to variable-disjoint 
factorization. 

Given a degree-d homogeneous polynomial / € F(A), we apply the following simple 
transformation to /: For each variable Xi € X we introduce d variables Xii,Xi 2 , ■ ■ ■ ,Xid. For 
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each monomial m € mon{f ), we replace the occurrence of variable Xi in the position of m 
by variable Xij. The new polynomial /' is in F({xjj}). The crucial property of homogeneous 
polynomials we use is that for any factorization f = gh both g and h must be homogeneous. 

Lemma 3.3 Let f € F(X) be a homogeneous degree d polynomial and f he the polynomial 
in¥{{xij}) obtained as above. Then 

• The polynomial f is variable-disjoint irreducible iff f is irreducible. 

• If f = 9 i 92 ■■■ 9t variable-disjoint factorization of f, where each g'^. is variable- 

disjoint irreducible then, correspondingly f = gig 2 ■ ■ ■ gt is a factorization of f into 
irreducihles g^, where g^ is obtained from g'^ by replacing each variable Xij in g'/^ by Xi. 

Proof. The first part follows because if / is reducible and f = gh then f = g'h', where g' is 
obtained from g by replacing the variables Xi by Xij, and h' is obtained from h by replacing 
the variables Xi by where s = deg 5 . 

For the second part, consider the product gig 2 ■ ■ ■ 9t- As all the factors gj- are homogeneous, 
it follows that each g^ is irreducible for otherwise g^ is not variable disjoint irreducible. 
Furthermore, any monomial m in gig 2 ... gt can be uniquely expressed as m = mim 2 ■ ■ ■ mt, 
where G mon{gk) for each k. Thus, for each i and j, replacing the occurrence of Xi by 
Xij in the product gig 2 ■ ■ ■ gt will give us f again. Hence gig 2 ... gt is f. □ 

It follows easily that factorization of homogeneous polynomials is reducible to variable- 
disjoint factorization and we can solve it efficiently using Theorems 12.71 and 12.81 depending 
on how the polynomial is input. We summarize this formally. 

Theorem 3.4 Homogeneous polynomials f G F(X) have unique factorizations into irre¬ 
ducible polynomials. Moreover, this factorization can be efficiently computed: 

• Computing the factorization of a homogeneous polynomial f given by an arithmetic cir¬ 
cuit of polynomial degree is polynomial-time reducible to computing the variable-disjoint 
factorization of a polynomial given by an arithmetic circuit. 

• Factorization of f given by an ABP is constant-depth reducible to variable-disjoint fac¬ 
torization of polynomials given by ABPs. 

• Factorization of f given in sparse representation is constant-depth reducible to variable- 
disjoint factorization of polynomials given by sparse representation. 


4 A Polynomial Decomposition Problem 

Given a degree d homogeneous noncommutative polynomial / G F(X), a number k in unary 
as input we consider the following decomposition problem, denoted by SOP (for sum of 
products decomposition): 

Does / admit a decomposition of the form 

/ = 51^1 H- Vgkhkf. 

where each gi G F(X) is a homogeneous polynomial of degree di and each hi G F(X) is 
a homogeneous polynomial of degree d 2 . Notice that this problem is a generalization of 
homogeneous polynomial factorization. Indeed, homogeneous factorization is simply the case 
when k = 1. 
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Remark 4.1 As mentioned in \Arl4- ^, it is interesting to note that for commutative polyno¬ 
mials the complexity of SOP is open even in the case k = 2. However, when f is of constant 
degree then it can he solved efficiently by applying a very general algorithm jArlf^ based on a 
regularity lemma for polynomials. 


When the input polynomial / is given by an arithmetic circuit, we show that SOP is in 
MA n coNP. On the other hand, when / is given by an algebraic branching program then 
SOP can be solved in deterministic polynomial time by some well-known techniques (we can 
even compute ABPs for the gi and hi for the minimum k). 

Theorem 4.2 Suppose a degree d homogeneous noncommutative polynomial f € F(^), and 
positive integer k encoded in unary are the input to SOP; 

(a) If f is given by a polynomial degree arithmetic circuit then SOP is in MAPlcoNP. 

(b) If f is given by an algebraic branching program then SOP is in deterministic polynomial 
time (even in randomized NC^J. 

(c) If f is given in the sparse representation then SOP is equivalent to the problem of 
checking if the rank of a given matrix is at most k. In particular, if¥ is the field of 
rationals, SOP is complete for the complexity class C=Lo 

We hrst focus on proving part (a) of the theorem. If {f,k) is a “yes” instance to SOP, 
then we claim that there exist small arithmetic circuits for the polynomials gi,hi, i G [k]. 

We define the partial derivative matrix Af for the polynomial / as follows. The rows of 
Af are indexed by degree di monomials and the columns of Aj by degree ^2 monomials (over 
variables in A). For the row labeled m and column labeled m', the entry is dehned as 

Am,m' = f{mm'). 

The key to analyzing the decomposition of / is the rank of the matrix Af. 

Claim 4.3 Let f G F(A) he a homogeneous degree d polynomial. 

(a) Then f can he decomposed as f = gihi + • • • -|- gkhk for homogeneous degree di poly¬ 
nomials gi and homogeneous degree d 2 polynomials hi if and only if the rank of Af is 
bounded by k. 

(b) Furthermore, if f is computed by a noncommutative arithmetic circuit C then if the 
rank of Af is bounded by k there exist polynomials gi,hi G F(A), i G [k], such that 
f = gihi + gkhk, where gi and hi have noncommutative arithmetic circuits of size 
poly{\C\,n,k) satisfying the above conditions. 

Proof of Claim \4.3^ 

For a homogeneous degree di polynomial g let g denote its coefficient column vector 
whose rows are indexed by degree di monomials exactly as the rows of Af. Similarly, for a 
homogeneous degree d 2 polynomial h there is a coefficient row vector [h with columns indexed 
by degree d 2 monomials (as the columns of Af). 

^The logspace counting class C=L captures the complexity of matrix rank over rationals |AB099 |. 
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Observe that if / can be decomposed as the sum of products f = gihi + ■ ■ ■ + Qkhk, where 
each Qi is degree di homogeneous and each hi is degree d 2 homogeneous then the matrix Af 
can be decomposed into a sum of k rank-one matrices: 


= 51^1^ H- 'rQkhk^. 

It follows that the rank of Af is bounded by k. Conversely, if the rank of Af is k then Af 
can be written as the sum of k rank 1 matrices. Since each rank one matrix is of the form gJi^, 
we obtain an expression as above which implies the decomposition of / as gihi gkhk- 

Now, if the rank of Af is k then there are degree di monomials mi,... ,mk and degree 
d 2 monomials m'^,... ,m^ such that the k x k minor of Af corresponding to these rows and 
columns is an invertible matrix K. W.l.o.g, we can write the p x q matrix Af as 


f K A2\ 

V Ai j; 

for suitable matrices Ai, A 2 , J. Moreover, since is rank k, we can row-reduce to obtain 


/ I 0\(KA2\^(KA2\ 

V -AiK-1 I J\A^ J J \0 O J 

and column reduce to obtain 

/ / 0\( K A2\( I -K-^A2 \^(K 0\ 

V -Aii^-i lylVAiJylVo / J \0 O J 

It is easy to verify that this yields the following factorization for Af = UlkV 


= 


K O 
Ai / 


I 0\( I K-^A2 \ 
O O )\0 I J 


Since we can write Ik as the sum eie'i eke'^ for standard basis vectors of suitable 

dimensions, we can express Af as the sum of k rank-one matrices (Uei){e'^V). We observe 
that the column vector Uci is the column of matrix U, which is identical to the 
column of matrix Af. Therefore, the polynomial represented by this vector corresponds 
to the partial derivative of / w.r.t the monomial m,, which can computed by a circuit of 
desired size. Similarly, any row vector e'^V is the row of matrix V, which is identical to 
the row of matrix Af, scaled by the transformation K~^. Moreover, the entries of K~^ 
are efficiently computable, and have size bounded by polynomial in the input size for the 
fields of our interest. Therefore, the polynomial represented by this vector corresponds to a 
linear combination of the partial derivatives of / w.r.t the monomials m [,..., m'^, which can 
computed by a circuit of desired size. Therefore, there exist polynomials gi,hi for i G [A:] 
which satisfy the conditions of the second part of the claim. □ 


Proof of Theorem 1^.2. 


Part (a). We first show that SOP is in MA. Given as input polynomial / G F(A) by 
an arithmetic circuit, the MA protocol works as follows. Merlin sends the description of 
arithmetic circuits for the polynomials gi, hi,l < i < k. By the second part of Claim this 
message is polynomial sized as the polynomials gi and gi have circuits of size poly{\C\,n,k). 


16 



Arthur verifies that / = gihi + 52^2 + • • • + using a randomized noncommutative PIT 
algorithm with suitable success probability. This establishes the MA upper bound. 

Next, we show that the complement problem SOP is in NP. I.e. given as input a polyno¬ 
mial / G F(A) such that / is not in SOP, we show that there is a poly{\C\,n, fe)-sized proof, 
verifiable in polynomial time, that / cannot be written as f = gihi + ■■■+ gi^hk- It suffices to 
exhibit a short proof of the fact that the rank of matrix Af is at least k + 1. This can be done 
by listing k + 1 degree di monomials mi ,..., mk+i and degree ^2 monomials m'l ,..., 
such that the {k + 1) x (A; -I- 1) minor ol Af indexed by these monomials has full rank. This 
condition can be checked in polynomial time by verifying the determinant of the minor to be 
non-zero. 

Part (b). When the polynomial is given as an ABP P, we sketch the simple polynomial-time 
algorithm for SOP. 

Our goal is to compute a decomposition 

/ = 51^1 + ■■■ + gkhk (3) 

for minimum k, where deg(( 7 j) = di and deg{hi) = d 2 - In order to compute this decom¬ 
position, we can suitably adapt the multiplicity automaton learning algorithm of Beimel et 
al |BB-|-00| . We can consider the input homogeneous degree-d polynomial / € F(A) as a 
function / : X'^ —>■ F, where f{m) is the coefficient of m in the polynomial /. The learn¬ 
ing algorithm works in Angluin’s model. More precisely, when given black-box access to the 
function f{m) for monomial queries m, it uses equivalence queries with counterexamples and 
learns the minimum size ABP computing / in polynomial time. Given a hypothesis ABP 
P' for /, we can simulate the equivalence queries and finding a counterexample by using the 
Raz-Shpilka PIT algorithm [RSOb] on P — P' . The minimized ABP that is finally output by 
the learning algorithm will have the optimal number of nodes at each layer. In particular, 
layer di will have the minimum number of nodes k which will give the decomposition in 
Equation [3l Furthermore, the ABPs for the polynomials gi and hi can also be immediately 
obtained from the minimum size ABP. 

Part (c). When the polynomial / G F(A) is given in sparse representation, we can explicitly 
write down the partial derivative matrix Af and check its rank. Moreover, we can even 
compute the decomposition by computing a rank-one decomposition of the partial derivative 
matrix. The equivalence arises from the fact that the rank of a matrix A is equal to the 
minimum number of summands in the decomposition of the noncommutative polynomial 
ttijXiXj as a sum of products of homogeneous linear forms. 

ije[n] 

□ 


An NP-hard decomposition problem 

We now briefly discuss a generalization of SOP. Given a polynomial / G F(X) as input 
along with k in unary, can we decompose it as a fc-sum of products of three homogeneous 
polynomials: 

/ = aiftici -I- a2b2C2 -|--h akbkCk, 

where each Oj is degree di, each bi is degree ^ 2 , and each c* is degree ds? 
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It turns out that even in the simplest case when / is a cubic polynomial and the ai,bi, Ci 
are all homogeneous linear forms, this problem is NP-hard. The tensor rank problem: given 
a 3-dimensional tensor Aijk, 1 < < n checking if the tensor rank of A is bounded by 

k, which is known to be NP-hard |Has90j is easily shown to be polynomial-time reducible to 
this decomposition problem. 

Indeed, we can encode a three-dimensional tensor Aij^ as a homogeneous cubic noncom- 
mutative polynomial / = Yli j A:G[n] ^ijkXiVjZk, such that any summand in the decomposition, 
which is product of three homogeneous linear forms, corresponds to a rank-one tensor. This 
allows us to test whether a tensor can be decomposed into at most k rank-one tensors, which 
is equivalent to testing whether the rank of the tensor is at most k. 

5 Concluding Remarks and Open Problems 

The main open problem is the complexity of noncommutative polynomial factorization in the 
general case. Even when the input polynomial / G F(X) is given in sparse representation 
we do not have an efficient algorithm nor any nontrivial complexity-theoretic upper bound. 
Although polynomials in F {X) do not have unique factorization, there is interesting structure 
to the factorizations [Co85[ IC^ which can perhaps be exploited to obtain efficient algorithms. 

In the case of irreducibility testing of polynomials in F(X) we have the following obser¬ 
vation that contrasts it with commutative polynomials. Let F be a fixed hnite field. We note 
that checking if / G F(A) given in sparse representation is irreducible is in coNP. To see this, 
suppose / is s-sparse of degree D. If / is reducible and f = gh is any factorization then 
each monomial in 5 or h is either a prefix or a suffix of some monomial of /. Hence, both g 
and h are sD-sparse polynomials. An NP machine can guess g and h (since coefficients are 
constant-sized) and we can verify if f = gh in deterministic polynomial time. 

On the other hand, it is an interesting contrast to note that given an s-sparse polynomial / 
in the commutative ring F[xi, X 2 , ■ ■ ■, we do not know if checking irreducibility is in coNP. 
However, checking irreducibility is known to be in RP (randomized polynomial time with one 
sided-error) as a consequence of the Hilbert irreducibility criterion [Ka89] . If the polynomial 
/ is irreducible, then if we assign random values from a suitably large extension field of F 
to variables X 2 ,... ,Xn (say, xi ^ ri) the resulting univariate polynomial f{xi,r 2 , ■ ■ ■, rn) is 
irreducible with high probability. 

Another interesting open problem that seems closely related to noncommutative sparse 
polynomial factorization is the problem of finite language factorization [SYOO] . Given as 
input a finite list of words L = {wi,W 2 ,---,Ws} over the alphabet X the problem is to 
check if we can factorize L as L = L 1 L 2 , where Li and L 2 are finite sets of words over X 
and L 1 L 2 consists of all strings uv for u G Li and u G ^ 2 - This problem can be seen as 
noncommutative sparse polynomial factorization problem where the coefficients come from 
the Boolean ring {0,1}. No efficient algorithm is known for this problem in general, neither 
is any nontrivial complexity bound known for it [SYOO] . On the other hand, analogous to 
factorization in F(X), we can solve it efficiently when L is homogeneous (i.e. all words in L 
are of the same length). Factorizing L as L 1 L 2 , where Li and L 2 are variable-disjoint can 
also be efficiently done by adapting our approach from Section [2j It would be interesting if 
we can relate language factorization to sparse polynomial factorization in F {X) for a field F. 
Is one efficiently reducible to the other? 
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